Protecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations

Molecular dynamics (MD) is one of the most common tools in computational chemistry. Recently, our group has employed accelerated molecular dynamics (aMD) to improve the conformational sampling over conventional molecular dynamics techniques. In the original aMD implementation, sampling is greatly improved by raising energy wells below a predefined energy level. Recently, our group presented an alternative aMD implementation where simulations are accelerated by lowering energy barriers of the potential energy surface. When coupled with thermodynamic integration simulations, this implementation showed very promising results. However, when applied to large systems, such as proteins, the simulation tends to be biased to high energy regions of the potential landscape. The reason for this behavior lies in the boost equation used since the highest energy barriers are dramatically more affected than the lower ones. To address this issue, in this work, we present a new boost equation that prevents oversampling of unfavorable high energy conformational states. The new boost potential provides not only better recovery of statistics throughout the simulation but also enhanced sampling of statistically relevant regions in explicit solvent MD simulations

On the Use of Accelerated Molecular Dynamics to Enhance Configurational Sampling in Ab Initio Simulations

We have implemented the accelerated molecular dynamics approach (Hamelberg, D.; Mongan, J.; McCammon, J. A. J. Chem. Phys. 2004, 120 (24), 11919) in the framework of ab initio MD (AIMD). Using three simple examples, we demonstrate that accelerated AIMD (A-AIMD) can be used to accelerate solvent relaxation in AIMD simulations and facilitate the detection of reaction coordinates: (i) We show, for one cyclohexane molecule in the gas phase, that the method can be used to accelerate the rate of the chair-to-chair interconversion by a factor of ∼1 × 105, while allowing for the reconstruction of the correct canonical distribution of low-energy states; (ii) We then show, for a water box of 64 H2O molecules, that A-AIMD can also be used in the condensed phase to accelerate the sampling of water conformations, without affecting the structural properties of the solvent; and (iii) The method is then used to compute the potential of mean force (PMF) for the dissociation of Na−Cl in water, accelerating the convergence by a factor of ∼3−4 compared to conventional AIMD simulations.(2) These results suggest that A-AIMD is a useful addition to existing methods for enhanced conformational and phase-space sampling in solution. While the method does not make the use of collective variables superfluous, it also does not require the user to define a set of collective variables that can capture all the low-energy minima on the potential energy surface. This property may prove very useful when dealing with highly complex multidimensional systems that require a quantum mechanical treatment

Ras Conformational Switching: Simulating Nucleotide-Dependent Conformational Transitions with Accelerated Molecular Dynamics

Ras mediates signaling pathways controlling cell proliferation and development by cycling between GTP- and GDP-bound active and inactive conformational states. Understanding the complete reaction path of this conformational change and its intermediary structures is critical to understanding Ras signaling. We characterize nucleotide-dependent conformational transition using multiple-barrier-crossing accelerated molecular dynamics (aMD) simulations. These transitions, achieved for the first time for wild-type Ras, are impossible to observe with classical molecular dynamics (cMD) simulations due to the large energetic barrier between end states. Mapping the reaction path onto a conformer plot describing the distribution of the crystallographic structures enabled identification of highly populated intermediate structures. These structures have unique switch orientations (residues 25–40 and 57–75) intermediate between GTP and GDP states, or distinct loop3 (46–49), loop7 (105–110), and α5 C-terminus (159–166) conformations distal from the nucleotide-binding site. In addition, these barrier-crossing trajectories predict novel nucleotide-dependent correlated motions, including correlations of α2 (residues 66–74) with α3-loop7 (93–110), loop2 (26–37) with loop10 (145–151), and loop3 (46–49) with α5 (152–167). The interconversion between newly identified Ras conformations revealed by this study advances our mechanistic understanding of Ras function. In addition, the pattern of correlated motions provides new evidence for a dynamic linkage between the nucleotide-binding site and the membrane interacting C-terminus critical for the signaling function of Ras. Furthermore, normal mode analysis indicates that the dominant collective motion that occurs during nucleotide-dependent conformational exchange, and captured in aMD (but absent in cMD) simulations, is a low-frequency motion intrinsic to the structure

Enhanced Lipid Diffusion and Mixing in Accelerated Molecular Dynamics

Accelerated molecular dynamics (aMD) is an enhanced sampling technique that expedites conformational space sampling by reducing the barriers separating various low-energy states of a system. Here, we present the first application of the aMD method on lipid membranes. Altogether, ∼1.5 μs simulations were performed on three systems: a pure POPC bilayer, a pure DMPC bilayer, and a mixed POPC:DMPC bilayer. Overall, the aMD simulations are found to produce significant speedup in trans–gauche isomerization and lipid lateral diffusion versus those in conventional MD (cMD) simulations. Further comparison of a 70-ns aMD run and a 300-ns cMD run of the mixed POPC:DMPC bilayer shows that the two simulations yield similar lipid mixing behaviors, with aMD generating a 2–3-fold speedup compared to cMD. Our results demonstrate that the aMD method is an efficient approach for the study of bilayer structural and dynamic properties. On the basis of simulations of the three bilayer systems, we also discuss the impact of aMD parameters on various lipid properties, which can be used as a guideline for future aMD simulations of membrane systems

Large-Scale Conformational Changes of Trypanosoma cruzi Proline Racemase Predicted by Accelerated Molecular Dynamics Simulation

Chagas' disease, caused by the protozoan parasite Trypanosoma cruzi (T. cruzi), is a life-threatening illness affecting 11–18 million people. Currently available treatments are limited, with unacceptable efficacy and safety profiles. Recent studies have revealed an essential T. cruzi proline racemase enzyme (TcPR) as an attractive candidate for improved chemotherapeutic intervention. Conformational changes associated with substrate binding to TcPR are believed to expose critical residues that elicit a host mitogenic B-cell response, a process contributing to parasite persistence and immune system evasion. Characterization of the conformational states of TcPR requires access to long-time-scale motions that are currently inaccessible by standard molecular dynamics simulations. Here we describe advanced accelerated molecular dynamics that extend the effective simulation time and capture large-scale motions of functional relevance. Conservation and fragment mapping analyses identified potential conformational epitopes located in the vicinity of newly identified transient binding pockets. The newly identified open TcPR conformations revealed by this study along with knowledge of the closed to open interconversion mechanism advances our understanding of TcPR function. The results and the strategy adopted in this work constitute an important step toward the rationalization of the molecular basis behind the mitogenic B-cell response of TcPR and provide new insights for future structure-based drug discovery

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics

An extension of the constant pH method originally implemented by Mongan et al. (J. Comput. Chem.2004, 25, 2038−2048) is proposed in this study. This adapted version of the method couples the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics, in an attempt to overcome the sampling issues encountered with current standard constant pH molecular dynamics methods. Although good results were reported by Mongan et al. on application of the standard method to the hen egg-white lysozyme (HEWL) system, residues which possess strong interactions with neighboring groups tend to converge slowly, resulting in the reported inconsistencies for predicted pKa values, as highlighted by the authors. The application of the coupled method described in this study to the HEWL system displays improvements over the standard version of the method, with the improved sampling leading to faster convergence and producing pKa values in closer agreement to those obtained experimentally for the more slowly converging residues

Using Selectively Applied Accelerated Molecular Dynamics to Enhance Free Energy Calculations

Accelerated molecular dynamics (aMD) has been shown to enhance conformational space sampling relative to classical molecular dynamics; however, the exponential reweighting of aMD trajectories, which is necessary for the calculation of free energies relating to the classical system, is oftentimes problematic, especially for systems larger than small poly peptides. Here, we propose a method of accelerating only the degrees of freedom most pertinent to sampling, thereby reducing the total acceleration added to the system and improving the convergence of calculated ensemble averages, which we term selective aMD. Its application is highlighted in two biomolecular cases. First, the model system alanine dipeptide is simulated with classical MD, all-dihedral aMD, and selective aMD, and these results are compared to the infinite sampling limit as calculated with metadynamics. We show that both forms of aMD enhance the convergence of the underlying free energy landscape by 5-fold relative to classical MD; however, selective aMD can produce improved statistics over all-dihedral aMD due to the improved reweighting. Then we focus on the pharmaceutically relevant case of computing the free energy of the decoupling of oseltamivir in the active site of neuraminidase. Results show that selective aMD greatly reduces the cost of this alchemical free energy transformation, whereas all-dihedral aMD produces unreliable free energy estimates

Infrared Multiple Photon Dissociation Action Spectroscopy and Theoretical Studies of Diethyl Phosphate Complexes: Effects of Protonation and Sodium Cationization on Structure

The gas-phase structures of deprotonated, protonated, and sodium-cationized complexes of diethyl phosphate (DEP) including [DEP − H]−, [DEP + H]+, [DEP + Na]+, and [DEP − H + 2Na]+ are examined via infrared multiple photon dissociation (IRMPD) action spectroscopy using tunable IR radiation generated by a free electron laser, a Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS) with an electrospray ionization (ESI) source, and theoretical electronic structure calculations. Measured IRMPD spectra are compared to linear IR spectra calculated at the B3LYP/6-31G(d,p) level of theory to identify the structures accessed in the experimental studies. For comparison, theoretical studies of neutral complexes are also performed. These experiments and calculations suggest that specific geometric changes occur upon the binding of protons and/or sodium cations, including changes correlating to nucleic acid backbone geometry, specifically P–O bond lengths and ∠OPO bond angles. Information from these observations may be used to gain insight into the structures of more complex systems, such as nucleotides and solvated nucleic acids

Allosteric effects in cyclophilin mutants may be explained by changes in nano-microsecond time scale motions

The relationship between molecular motion and catalysis in enzymes is debated. Here, simulations of cyclophilin A and three catalytically-impaired mutants reveal a nanosecond-scale interconversion between active and inactive conformations, orders of magnitude faster than previously suggested